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Abstract 

We propose a simple density functional expression for the upper bound of the kinetic energy 
for electronic systems. Such a functional is valid in the limit of slowly varying density, its validity 
outside this regime is discussed by making a comparison with upper bounds obtained in previous 
work. The advantages of the functional proposed for applications to realistic systems is briefly 
discussed. 
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I. INTRODUCTION 



Kinetic energy functional of the electron density have been since long a subject of intense 
investigation. Startingwith the pioneering work of March [1] through the monumental work 
of Lieb (see e.g. the aim was to construct functionals which are accurate enough 

to properly describe physical and chemical properties and at the same time simple enough 
to allow a feasible computation. There has been a period of intense activity around this 
subject across the seventies until the end of the eighties during which many interesting 
results where produced. After about a decade of weaker activity, the interest in the subject 
has got a new vigor (see e.g. the interesting work of Ludena and coworkers fj, Q, a E| 
and the topical review of March ^(|). The reason of this renewed interest lies in the fact 
that in the meanwhile novel computational schemes for quantum calculations, where the 
kinetic functional plays a key role, have been developed. Of particular interest is the linear 
scaling real space kinetic energy functional method, where the kinetic energy is calculated 
as a functional of the electron density. The electron wavefunctions are no more required and 
for this reason the method is called orbital-free density-functional theory (OFDFT) (see 
13, [l4, 15, JjJ). Since neither the diagonalization of the electronic Hamiltonian 



e.g. 
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nor the reciprocal space sampling are required, such techniques allows for studies of relatively 
large systems compared to those treatable with the standard Kohn-Sham based approaches. 
Moreover, the fact that the calculation are done only in real space allows for the development 
of efficient quantum-classical interfaces which in current research are highly desirable within 
the emerging multiscale modeling techniques. For this reason we turned the attention to the 
derivation of a simple and physically well founded kinetic functional. We start from the most 
general polar form of the electron wavefunction and derive an upper bound which is exact 
in the limit of slowly varying density; we discuss its validity beyond such an approximation 
by comparing our results with those available in literature. A first interesting result is 
that our functional is potentially a better bound compared to some of those available in 
literature and by now well established; it is also universal, i.e. does not show explicit 
dependence on N, the number of electrons, differently from most of those found in literature. 
Next we combine our upper bound with the well known lower bound of Lieb and Thirring. 
By doing so, we conclude that a valid functional is likely to have the form of a Thomas- 
Fermi- Weizsacker type (with different constants), where the multiplicative constant of the 
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Weizsacker term is the only free parameter. To our knowledge this is the first time that a 
universal functional containing a Thomas- Fermi-like term and the Weizsacker one is obtained 
within the same derivation, and not by separated schemes. Moreover, the link we find 
between the wavefunction phase factor and the kinetic functional suggests an operative way 
to estimate the unknown constant. 

II. THE KINETIC ENERGY 

Let us consider a system of N electrons in a 3iV-dimensional volume Q N . The general N- 
particle wavefunction in polar form is: ?/>(ri,r 2 , ....rjy) = #(ri,r 2 , ....rA^e^ 1 " 1 ' 1 " 2 '"" 1 "^* 1 , where 
9(ti, r 2 , ....tat) and S(ri, r 2 , ....tat) are real functions in Q N . We also require 9(ti, r 2 , ....tn) 
to be antisymmetric with respect to any pair permutation, i.e. 9(ri, r 2 , ..rj, ..r 3 -..rjv) = 
—9(rx,r 2 ,..rj,..r i ..r^) and S(ri, r 2 , ••■■rjv) symmetric. In this way the wavefunction 
^(ri, r 2 , ....rjy) is antisymmetric as should be for a system of fermions. For simplicity the 
spin variables are not explicitly considered. We define the one electron density p(r) in Q 
as: p(r) = iV J nJV _i [9(ti, ...r^_ 1; r i; r i+1 , ..r N )} 2 dri..dr i _ 1 dr i+1 ..dr N , where the index i is ar- 
bitrary and can take any value from 1 to N; for the case i — 1, dr^i is not considered while 
for i = N, dr i+ i is not considered. The one electron density, in turn, satisfies the condition: 
j n p(r)dr = N. We will use atomic units H — 1, the electron mass m = 1 and the electron 
charge e = 1. Let us consider the average kinetic energy for the state ip: 

= -\ [ fV^/w (1) 
^ Jn N 

where d N uj = 11^ dr, and V = Vj. The expression above can be written also as (see 
e.g. 0): 

T* = \l IW| W (2) 
Substituting the expression ip = 6e tS into Eq|2]one obtains: 

T^= / + ^T~ d u - ( 3 ) 

Jn N V 1 1 

The integrand on the r.h.s. of Eq|3]is a 3iV-dimensional function; the goal is to reduce it, 
as rigorously as possible, to a simple three dimensional functional of p(r). 
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A. Upper bound to \VS\ in the slowly varying density limit 



The average over f2 n of the 3iV-dimensional momentum vector of the system is defined 

as: 



(V>MV>|) n w = - / Im[i)*Vi)]d N uj (4) 
Jn N 

substituting the polar form of ip into Eq0] we obtain: 

f Im[ij*Vij}d N u = [ [i>*VSi/;}d N uj. (5) 

Through EqEJwe find the relation (■?/>|m|'?/>) QJ v = {i/j\VS\i/j) qN , which suggests an interpre- 
tation of VS as a 3iV-dimensional quantum velocity field. Such an interpretation would not 

) , above all in the context of 
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201 ]. If we interpret VS" as a quantum 



be new, and can be often found in literature (see e.g. 
a fluid dynamics formulation of quantum mechanics 
velocity field than, in the limit of p(r) being a slowly varying function, we can state the 
following inequality for the absolute value of V5: 

(i>\\VS\\i>)<(i>\\P F M) (6) 

which explicitly reads: 

„ JV N 

/ e 2 \Tvs\d N uj< / e 2 \Tp F .\d N uj. (7) 

Where P Fi = Pf(rj), with = YliLi P-F( r «)> is the Fermi momentum, i.e. the maximum 
momentum a particle can reach in the limit of slowly varying density. Since Eq|7] must hold 
for any arbitrary subvolume of fl N we obtain: 

|V5| < iP^hVRG^ (8) 

and thus 

|V5| 2 < |P FjV | 2 ;VRg Q N . (9) 

From the inequality El we have: 

/ (P\VS\ 2 i(>)d N u< f (r\^F N \ 2 ^)d N u (10) 
Jn N Jn N 



or equivalently: 



/ 9 2 \VS\ 2 d N u< / p(r)|P F (r)| 2 rfr. (11) 
Jq n Jn 



The second term on the r.h.s. of Eqllllis obtained in the following way: 

/ {p\p F \ 2 4)d»u, = f: [ |P Fi | 2 rfr 4 / ^-^ = yf n^l p F,l 2 *i- (12) 

Ju N i^[Jni Jn"- 1 i^i^i 7V 

If we go back to EqEl since: |Pf(i")| = Cp[p(r)]^ 3 , where Cf is a constant, (CV) 2 = (37r 2 ) 2 / 3 , 
we have: 

/ e^d»u+[ if>C^< [ P (r)? /3 + [ |f>0| 2 d 
Jn N 1 Jn N 1 ~i Jn 1 Jn N 1 ~[ 



to (13) 



In this way we have determined an upper bound to T^, where the first term is written as 
a functional of p(r). We shall now reduce also the other term into a three dimensional 
expression, possibly as a functional of p(r). One may wonder why to study the phase 
factor S for atoms and molecules in their ground state. Indirectly this work suggests an 
unconventional answer; strictly speaking S = for the ground state of atoms and molecules, 
and this would mean that, in principle, the only kinetic functional term is the Weizsacker 
term, as it is discussed by Sears et al. (2^, Herring j3] and Luo 

0- However in Refs.Q, 



is discussed and shown that this results would not be correct and there must be an 
additional term to the kinetic energy coming from some angular part of the wavefunction 
which factorizes, i.e. something related to a phase factor. At the same time if tp = 9e lS 
is a solution, also ip* = 9e~ lS is a solution with the same energy; this also holds for linear 
combinations of the two. For atoms and molecules one can always choose, without loss of 
generality, linear combinations of ip and ip* that lead to real wavefunctions. In this case, as 
it can be verified by a straightforward calculation, the r.h.s. of EqEDis still an upper bound 
to the kinetic energy. 

B. Weizsacker and non local Information functional 

The term f nN \\79\ 2 d N r = YliLi fn N \^i9\ 2 d N T is what Sears et al. (21] refer to as the 
multivariate kinetic functional and can be written as: 



^ttJn^ 8Jn P(r) 8 J QN N f(r 1 ,...r i _ 1 ,r i+1 ...r N /r i ) 

(14) 

where the first term on the r.h.s. is the well known Weizsacker term and i G [l,iV] with 
...r i _i,T i+ i...r^/r i ) is proportional to the conditional density, i.e. the probability den- 



sity of finding a certain spatial configuration for the N — 1 particles once the position of the 



i-th particle is assigned (see also 24j). Following the work of Ref.j2l| we can write: 

N 

2 



where 



/ p(r)I(r)dr = / p(r) 



\Vif{r 2 ,r 3 ....r N /r n 

-dr2dr 3 ...dr^ 



2 

' dr (16) 



-i /(r 2 ,r 3 ,...rjv/r) 

with J(r) being the well known non local information functional within Fisher information 
theory. An exact expression for 7(r) is difficult to find, however one can notice that I(r) 
would be equivalent to write in local form some kinetic correlation effects which are usually 
considered negligible Q|. Combining the results above with those of the previous section we 
obtain: 

r ,^/ nWr)]5 /3 + i| n ^ r+ iy /(r)/(r)dr (17) 

III. BEYOND THE SLOWLY VARYING DENSITY REGIME 

In this section, by comparing the result of EqJTTl with upper bounds to available in 
literature, we will discuss the validity of our results for the general case. For simplicity we 
will neglect the non local information functional term since, as said before, it represents a 
non relevant correction. Under this hypothesis Eq lT7l becomes: 

It is encouraging to note that the same functional form for the upper bound, with a slightly 
different multiplicative constant for the first term, was conjectured, following arguments 
different from ours, by Lieb [4|. March and Young obtained a result similar to ours but the 
constant multiplying the J p 5 ^ 3 dr must be determined for each atom [lj]. Gazquez and Robles 
|25j derived a kinetic energy functional of the form: = d (l - A) / p 5 / 3 rfr+| / l -^dr, 
the same functional was independently proposed by Acharya et al. |26j . This functional has 
explicit dependence on N, the number of particle, thus it is not universal, however, in the 
thermodynamic limit, i.e. N large, has the same functional form of the functional we de- 
rived (with, again, a slightly different constant for the first term). Moreover our result is 



fully consistent with the estimate of the relative-phase- energy term of Herring 



221 where he 
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concludes that the upper bound must be something which has the same form as our func- 
tional. Even better comparison can be made for, the simpler, one dimensional case; now the 
term J p 5 ^ 3 dr in Eq J18I becomes J p 3 dx because in our derivation Pp ~ p in one dimension. 
The resulting upper limit has the same functional form for obtained by Harriman 27 1 
using the special equidensity orbitals (SEDOs) construction for ip (see also 28]). Moreover 
it represents an upper bound, in the limit of large N, to the rigorous functional found by 
March and Young in one dimension: < const x J p 3 / 2 dx + J ^-^-dx. The arguments 
above, although do not represent an explicit proof, suggest that the functional we propose 
may indeed be a valid upper bound beyond the slowly varying density approximation. An 
intuitive argument to support the validity of our hypothesis is the following. In the slowly 
varying density regime, our upper bound condition on the momentum of the single electron 
is rather " large" . In fact we approximate the momentum of each electron with the maximum 
value possible instead of distributing among the N electrons all the states available from 
|m| = to |m| = |Pjr|. By slowly moving from this regime to an intermediate one, for a 
large number of particles, despite the fact that Pp slowly looses its physical meaning, the 
upper bound hypothesis |Vi>| < const x p 1 / 3 it is likely to still hold for an extended range of 
densities. Moreover, the larger N the more extended the range of validity; in the thermody- 
namic limit one can expect such an upper limit to be always valid, as the comparison with 
the functional of Gazquez and Robles [25] and of Acharya et al. [26] suggests. Interestingly, 
if this was the case, we would have found a better upper bound to compared to that of 
Zumbach j^: 

T<p < [I + C Zu N^} 1 - J l -^dr; C Zu = 15(4tt) 2 ^ QV^ (19) 

which was obtained using a SEDOs-like construction of ip and is considered of general 
validity, although its explicit dependence on N makes it not universal. In fact, as it is 
discussed by Pathak and Gadre [30j] , and also reported by Zumbach [29] , by using Schwartz 
inequality and then applying one Sobolev inequality in three dimensions it is possible to 
obtain the following relation: 

J p 5 ' 3 dr < C PG x iV 2 / 3 i J ^!rf r (20) 
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where Cpg is a constant (see |30(), if now we compare EqJTHlwith EqUHlvia Eq|20J we obtain: 



[C_ 12 

since Cp G x Cp < C^ u . 



i/^^+^W/W* (21) 

% Jn P ° J P 



IV. UPPER AND LOWER BOUND FOR T, 



An interesting property of our upper bound can be obtained by relating it to the Lieb 



Thirring inequality and to the consequent lower bound to T, 



ea i iv r 



31. 



^ > % / P 5/S rfr. (22) 



From Refs. 



3J it follows that C£ of EqUHJ is very close to C LT of Eq|22] In fact 



C'p ~ 9.57 while Clt ~ 9.11. Actually Lieb in Ref.J^J argues that numerical calculations 
improve Clt to 9.578. If we consider Cp ~ Clt, then we can write in good approximation 
the following relation: 

C x J p^ z dr <T^<C x y p 5 / 3 rfr ^y~ dr ( 23 ) 

where C is a constant whose value is between that of Cp/2 and Clt/2, thus very close to 
both. The implications of Eq|2H]are very interesting; it suggests that a valid approximation 
for can be written as: 

T^Cx / p b ' z dv + qxl f ^^dr (24) 
J 8 Jq P 

where < q < 1, is the only free parameter. A possible way to obtain a first guess for q 

would be by estimating VS 1 numerically for some simple model systems. However already 

at this stage Eq|2H represents a very interesting result, as it has been previously discussed 

in the introduction, in connection with current quantum based computational schemes. 



V. DISCUSSION AND CONCLUSIONS 

We have used the polar form of a many-particle electron wavefunction to derive a simple 
functional for the kinetic energy. Such a functional represents an upper bound to the true 
kinetic energy and it is exact in the limit of slowly varying density. We have discussed its 
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validity beyond such an approximation by making a comparison with functionals obtained 
in previous work. The novelty of our derivation presents different aspects; we obtain within 
the same theoretical framework a kinetic functional containing both a Thomas-Fermi like 
functional and the Weizsacker functional. These are usually determined by following sep- 
arated procedures under different physical approximations. This result provides a physical 
justification to a high desirable upper bound functional which was heuristically conjectured 
by Lieb -4]. Another advantage is that it does not show explicit dependence on N , differ- 
ently from most of the upper bounds available in literature, and represents a better bound 
compared to other well established functionals. Finally, combined with the Lieb-Thirring 
lower bound leads to the conclusion that a valid functional, which can well approximate the 
true one, could be determined by simply tuning the multiplicative constant (between zero 
and one) of the Weizsacker term. This result is very interesting above all for applications 
to condensed matter systems within free orbital density functional schemes. In fact, it pro- 
vides not only a theoretical background to justify the kinetic functional currently employed, 
but also a possible procedure to determine new ones. Approaches as the one shown in this 
paper, are crucial for the development of computationally efficient and theoretical flexible 
quantum mechanical techniques for modern multiscale simulations . As discussed into the 
introduction, the computational apparatus is available and new ideas, which can improve 
the energy functionals currently available or suggest new ones, are strongly required. In this 
respect, our contribution suggests a possible way to proceed. 
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